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The OH molecule is currently of great interest from the perspective of ultracold chemistry, quan- 
tum fluids, precision measurement and quantum computation. Crucial to these applications are 
the slowing, guiding, confinement and state control of OH, using electric and magnetic fields. In 
this Letter, we show that the corresponding eight-dimensional effective ground state Stark-Zeeman 
Hamiltonian is exactly solvable. Based on our solution, we analyze a waveguide for OH molecules. 
Our analytical result opens the way to insightful characterization of the magneto-electrostatic ma- 
nipulation of ground state OH, and is also immediately applicable to other diatomic free radicals 
such as CIO, SF, NS, OD and SrF. 

PACS numbers: 33.20.-t, 33.15.Kr, 37.10.Pq 



The OH molecule in its ground X^Il3/2 state is 
presently widely employed in investigations of ultracold 
chemistry [l|-[3|, precision measurements and quan- 

tum computation Q]. Particularly interesting is the re- 
cently implemented evaporative cooling of OH close to 
Bose- Einstein condensation With such experiments 
underway, the exploration of quantum degeneracy and 
molecular optics [9[ with OH should shortly become re- 
ality. 

A substantial reason behind the suitability of OH as a 
workhorse for the mentioned experiments is the fact that 
it is a polar paramagnetic molecule, i.e. it carries both 
electric and magnetic dipole moments. Electric and mag- 
netic fields can therefore be used to slow, guide, confine 
and generally manipulate OH It follows that a 

quantitative as well as qualitative understanding of the 
corresponding Stark-Zeeman spectrum is of great rele- 
vance. 

In this Letter we present the exact solution of the 
eight-dimensional Stark-Zeeman Hamiltonian of OH in 
its ground state [l3|- This molecular Hamilto- 

nian is an effective one, neglects hyperfine structure, and, 
to the best of our knowledge, has only been solved numer- 
ically [i, 12, 1^. Exact analytical solutions to a physical 
problem are of value because they describe how the sys- 
tem behaves in any regime, as opposed to numerical solu- 
tions, which usually describe a specific parameter range. 
During the preparation of this article, an analytical solu- 
tion of the field-depende nt p art of the Hamiltonian was 
reported by Bohn et al. [iSj . 

To demonstrate the usefulness of our solutions, we an- 
alyze a simple magneto-electrostatic waveguide for OH, 



consisting of a current-carrying wire, a magnetic bias 
field, and an electric offset field [l6|. We calculate analyt- 
ically the energy minimum, trapping frequency, and an- 
harmonicity of the waveguide as a function of the electric 
field. Using realistic parameters [1, 12, 13 1, we demon- 
strate the advantage of using an electric field in place of 
the standard magnetic offset field. 

We expect that the solutions presented in this Let- 
ter can be used to construct an analytical yet realistic 
theory of nonadiabatic processes in traps, which so far 
has relied on the four-dimensional model of a fictitious 
^Hi/2 OH ground state 11 71. Our work may also be of 
relevance to atmospheric [18|, interstellar [19| and com- 
bustion physics [20|, where OH plays an important role; 
and is immediately applicable to other molecules which 
have ^H3/2 Hund's case (a) states in their spectra, such 
as CIO, SF, NS, OD and SrF. Lastly, we hope that our 
discovery of an exactly solvable model will be of inter- 
est to physicists investigating the behavior of atoms and 
molecules in strong fields, and whose concerns include 



integrability and quantum chaos [21 



We begin with the Stark-Zeeman Hamiltonian for OH 
in the state, which has been derived earlier 113 1 
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where Ho is the field-free Hamiltonian, and /ih are the 
electric and magnetic dipole moments of the molecule, 
respectively, and E [B] is the electric [magnetic] field im- 
posed on the molecule. It has been shown that this model 
is valid when hyperfine structure is negligible, such as for 
electric fields stronger than IkV/cm, and magnetic fields 



above lOOG For the present work, we will assume OH 
vapor temperatures higher than 2.7mK, in which case hy- 
perfine structure can also be neglected [s, 13|. 

The matrix representation of the Hamiltonian in 
Eq. ([1]) can be obtained using the Hund's case (a) parity 
basis I J, M, e) suggested by Lara et ai, where J — 3/2 
is the total angular momentum, M its projection in the 



laboratory frame, f2 its projection on the internuclear 
axis, and e = {e, /} is the parity 17 1. Following Ref . [TtI . 
both the electric and magnetic moments are assumed to 
lie along the axis of the molecule, and the magnetic field 
is chosen along the laboratory z axis, with which the elec- 
tric field makes an angle 9. With these assumptions, the 
Hamiltonian matrix has been found to be 13 1 
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where A is the lambda-doubling parameter, is the 
Bohr magneton, fj,e the molecular electric dipole moment, 
and E [B] are the magnitudes of the electric [magnetic] 
fields. 

To the best of our knowledge, the matrix Hm of 
E q.(Et has only been diagonalized numerically so far 
[sl-lli. 13, 1^. We now show, remarkably, that this eight- 
dimensional Hamiltonian can be solved analytically. This 
can be seen by calculating the characteristic polynomial 
P(A) of Hm, 



P{X) = \Hm- XI\ = J2Pny 



(3) 
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The coefficients of all the odd powers of A vanish: pi = 
P3 = P5 = P7 = 0. Thus, -P(A), which is an octic in A, 
can be written as a quartic in A^, 

P(A) = po+P2X^+Pi{XT+P6{X^)'+P8iX^)\ (4) 



and therefore its roots can be found analytically [22] ■ Ex- 
pressions for the coefficients Pn and eigenvalues i?M,£ of 
Eq. ^ arc most readily obtained using Mathematica and 
are provided in the Supplementary Material. The eigen- 
value of interest in this article, -E3/2,/, can be written 
as 
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We have verified that the analytical eigenvalues, dis- 
played in Fig. [1] reproduce exactly t he p urely numeri- 
cally obtained OH spectra, as in Ref.[13|, for example. 
We note that although it is conventional and useful to 
label the eigenvalues using the angular momentum pro- 
jection M and the parity e, neither is a good quantum 
number [l^. 

As expected from Eq. ^ and shown in Fig.[Tl the spec- 
trum is vertically symmetric about zero energy, i.e. if A is 
an eigenvalue, so is — A. The observation of this symmetry 
in the numerically generated spectra [13[ was the crucial 
step in arriving at the analytical solution of the problem. 
This kind of symmetry seems to have been noticed in 
only a handful of physical systems . Perhaps the best 
known example is that of a free Dirac particle, where 
the 'reflection' symmetry of the spectrum follows from 
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FIG. 1: Plot of the analytical eigenvalues of OH with the 
parameters A = 27r x 1.667GHz, fie = 1.66D, E = 2kV/cm, 
and 9 = n/2. The horizontal axis denotes the magnetic field 
in Tesla and the vertical axis corresponds to energy in units 
of temperature. For clarity, only the state E^^^j currently 
used for trapping OH [H, [3 has been labeled. 



the anticommutation of the charge conjugation operator 
with the corresponding Hamiltonian 2J] . In the present 
case, the physical origin of this spectral symmetry is not 
clear, partly due to the effective nature of the Hamilto- 
nian, and is an open and interesting question posed by 
our work. 

We now demonstrate how our solutions can be used 
to analytically characterize magneto-electrostatic ma- 
nipulation of OH molecules. The dramatic influence 
of electric fields on magnetic OH traps has recently 
been demonstrated experimentally, in terms of increasing 
steady-state trap population [l^l, as well as implement- 
ing evaporative cooling Q. These experimental configu- 
rations are quite complicated and will be considered in 
the context of our analytical solutions in future work. In 
the present article, we examine instead a simpler geom- 
etry. Particularly, we are interested in the effect of an 
electric field on a magnetic waveguide for OH. Such an 
electric field could be useful for polarizing the OH electric 
dipole, or for manipulating cold collisions [H, [Ij. 

We first briefly review the standard magnetic 
quadrupole waveguide configuration, which can be used 
to guide any paramagnetic molecule, including OH [l6| . 
As shown in Fig. [21 the guide consists of a wire carrying 
a current / in the z direction (i.e. out of the page), and 
a homogeneous magnetic bias field Bt in the x direction. 
We have chosen this configuration since it is analytically 
simple to model. The resulting magnetic field away from 
the origin is given by 



B, 



27r(a;2 + y2) 



y, (7) 



and has a zero at Xo ~ and yo = fioI/{2T:Bb). Para- 



FIG. 2: Schematic of the OH magneto-electrostatic guide con- 
sidered in this work. Included in the diagram are a wire carry- 
ing current I along the z axis out of the page, and a magnetic 
bias field Bi, along the x axis. Their superposition results in a 
quadrupole guide at the location marked {xo,yo)- Any point 
in the guide can be labeled using the coordinates p and (j), as 
shown. In the standard configuration, a magnetic offset field 
Bo along the z axis is also used. We consider in this work 
instead an electric field E, applied in the same direction. 



magnetic atoms and molecules can be guided along this 
line, parallel to the z axis. Shifting the origin to the 
point {Xo,yo), the magnitude of the magnetic field can 
be written as [id 



BbP 



yjp^ +yl + 2pyo cos(^' 



(8) 



where the radial distance p and the angle <f> are mea- 
sured as shown in Fig. [2] For p <C j/o, Eq. ([8]) is radially 
symmetric about {xo,yo)- Thus, we may simply write 
\Bt \ ~ bp, where b — Bh/yo is the magnetic field gradient 
at the guide center. 

One of the problems with this configuration is that 
due to the zero field at the trap center {p = 0) Ma- 
jorana losses can occur [l6l |. In order to prevent such 
losses, usually a magnetic offset field Bo along the z axis 
is applied, resulting in a total field magnitude B{p) = 
(|i?t|2 + B2)1/2. Assuming, as in current experiments, 
that the i?3/2,/ internal state is used for guiding OH [l3| . 
the potential energy of the molecule near the guide cen- 
ter can be written as a Taylor expansion of the exact 
eigenvalue 



E3/2AP) = 4%ji^) + Eiy,j{0)p + E',%j{0)p' + 
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In Eq. ®, 4;^,/0) = (n!)-ia"£;3/2,/(p)/ap"|p=o. Fur- 
thcr, the energy depends only on p since the system is 
cylindrically symmetric about, and translationally invari- 
ant along, the z axis. Lastly, except at the guide center, 
the net magnetic field does not point along the z axis, 
thus apparently contravening the convention of Eq. ([2|). 
However, in that case a new z axis can readily be cho- 
sen in the direction of the local magnetic field, and the 
procedure yields the same results as presented below. 

The required eigenvalue for this case can be found di- 
rectly by setting £' = in Eq. ([2]), which then takes a 
diagonal form. Using the resulting solution i?3/2,/(p) — 
+ 6pbB{p)/5, we find that -E^/2j(0) = hA/2 + 
6pbBo/5, implying that the nonzero magnetic field at the 
guide center interacts with the OH molecule, thus lower- 
ing Majorana losses. Further, we find £^3)2./(0) === and 



^3/2 f(^^ ^ SfJ-sb^ /5Bo, making the trap harmonic to 



3/2,/ 



lowest order in p. Lastly, Ey]^ /(O) =0 and Ey'^ y(0) = 
—8/156^/205^. Thus, the lowest order, i.e. quartic, an- 
harmonicity cannot be fully removed, although it can be 
reduced in magnitude by choosing a small b and a large 
Bo. 

We now consider the effect of replacing the magnetic 
offset field Bo by a homogeneous electric field E, also 
applied in the z direction 13 1, which conserves the 
cylindrical symmetry of the system. Using a realistic 
simulation, we have confirmed that for the parameters 
presented in this Letter the electric field near the guide 
center is only negligibly perturbed by the presence of the 
wire. The quantities relevant to guiding can be found 
from Eq. © with E 0,9 = 7r/2, and Eq. ([9]); such 
a treatment is equivalent to choosing a molecular z axis 
along the direction of the local magnetic field. 

We find that the energy at the guide center is given 
exactly by 

9 1/2 



(4) 



E. 
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The non-zero electric field at the trap minimum thus pro- 
vides the molecule with a quantization axis and may be 
expected to mitigate Majorana losses 15|. We note that 



the first derivative in Eq. (|9]) vanishes, i.e. £'3/2 /(O) = 0, 
making the trap harmonic to lowest order in p. The 
guide frequency = [2E^^\ ^{0)/m]^^^/2TT, where m is 
the mass of the OH molecule, is found to be 
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and is plotted in Fig. [H^a) for realistic parameters. 



(11) 




(b) 



5 10 15 20 
E (kV/cm) 



FIG. 3: (a) The frequency and (b) anharmonicity of the guide 
in Fig[2] as a function of the electric field. Parameters used 
are I = 1500A, Bb = 0.03T, m = 17amu, A = 27r x 1.66GHz, 
and He = 1.667D [l^ . implying yo = 10mm and b — 3T/m. 



Next, we examine the trap anharmonicity. The 
cubic anharmonicity is identically zero, while the 
quartic contribution ap can be found analytically 
to be ap = F/G, where F = 3(5;iA)'^ + 
57{5hA)^{2peE)^+216{5hA)^{2peE)-^~6A8{2^ieE)^, and 
G = ^{2neEf[9{2neE)^ + {5hAff/^. The denominator 
G is always positive, but the presence of the last nega- 
tive term in the numerator F immediately suggests that 
the anisotropy can be removed by choosing the electric 
field magnitude appropriately. Indeed, Fig. ^h) shows 
that the OH guide anharmonicity vanishes at about 
E = 3.5kV/cm. Choosing this value of the field might be 
useful as the effect of the anharmonic term can be quite 
large; at i? = 500V/cm, and p = 2mni it equals 20% that 
of the harmonic potential energy. Finally, for complete- 
ness, the depth of the guide as a function of the electric 
field has been provided in the Supplementary Material. 

Before concluding, we note two subtle points. First, 
at the center of the guide, B = 0, leading to 9 being 
undefined in Eq. However, this is not a problem, 

as Eq. ([TU| can be shown to hold for any 9, and the ex- 
pressions for the frequency and anharmonicity presented 
above can be recovered in the limit of arbitrarily small 
(but nonzero) B close to the guide center. Second, we 
note that both the guide frequency and anharmonicity 
diverge as E' — > 0. This is due to the avoided crossing 
turning into a crossing, and introducing a discontinuity 
in the derivatives of the Taylor expansion of Eq. ^ . At 
E ^ 0, the transverse character of the guide changes fun- 
damentally, from harmonic to linear. 

In conclusion, we have solved exactly the effec- 
tive Hamiltonian of the X2H3/2 state Hund's case (a) 
molecule in combined electric and magnetic fields, ne- 
glecting hyperfine structure. In the context of a realizable 
waveguide, we have shown that our results allow insight- 
ful analytic characterization of the magneto-electrostatic 
manipulation of OH molecules. Our analysis can read- 
ily be extended to other field configurations [0 , 12 1 and 
other molecules such as CIO, SF, NS, OD and SrF. 
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Below we provide additional details of the calculations presented in the main article. The analytical as well as 
numerical calculations were carried out using Mathematica; a notebook containing these results is available upon 
request to mxbsps@rit.edu. 

A. Coefficients of the ciiaracteristic polynomial 

The coefficients p„ in Eq.(3) of the main article arc found to be 

Pi = P3 =P5=P7 =0, (1) 



and 



^xmpf'K^ + 1185* A"* - 264:I3^E^A* + llsi^^A"* - 2013^ + 20E^A^ + A^ 

-AB^E^ (SIB^ + 81E'^ + 54^2 - 54E^A^ - 7A*) cos 29 + IQ2B^E^ cos 46l| , (2) 

P2 = (-95^ - 9^*^ - ^A^ - —E^K" - 3E^A* - 9B^E^ - —B*A^ - 9B^E^ 

^ 50000 \ 5 5 5 

+B^A* + ^B^E^A:^ + 2I3'^E^ ^95^ + 9i?2 + yA^ j cos26'l , (3) 



Pi 



— \59B^ + iQB^E^ + lOB^ A^ - 82B^E^ cos 261 + b9E^ + 30l;^ A^ + 3A''| , (4) 
5000 L i 



P, = -\\b'^ + E' + \K^Y (5) 
P8 = 1. (6) 

where B = AjisB, E = 2iJ.eE, and A = 5fi.A. (Continued on the next page). 
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B. Eigenvalues 

The roots of the characteristic polynomial of Eq. (4) in the main article are found most directly in analytic form 
using Mathematica, and are given below: 



EV2. = \/-^ + ^-^^^^, (8) 
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The Taylor series expansions provided in the article were obtained using Mathematica and the validity of the Taylor 
series up to the required number of terms was established by comparison with the full numerical solution for specific 
parameters. 



C. Reflection symmetry in the spectrum 

In the spectrum of the OH molecule, if A is an eigenvalue, so is —A, as shown in Fig.l in the main article. In many 
systems that exhibit such a feature, the symmetry of the spectrum follows from the existence of an operator C that 
anticommutes with the system Hamiltonian Hg, i.e. 

CHs + HsC = 0, (16) 

If we consider an eigenvector ip^ with an eigenvalue E, we can write Hs'4'+ ~ Multiplying from the left by C 

and rearranging, we find CHs4>+ = ECt]-J+. Finally, using the anticommutation, we arrive at 

Hs{Ccp+) = -E{Ccb+). (17) 

This implies that Ccjjjf. is an eigenfunction of Hs with an eigenvalue —E. The existence of the operator C thus leads 
to the ±E pairing of eigenvalues in the energy spectrum. Assuming C is invertible we can write Eq. (|16p as 



CHsC-^ = -Hs, 



(18) 



3 



25 
20 

s 

Q° 10 
5 


2 4 6 8 10 

E ( kV/cm ) 

FIG. 1: The transverse depth of the guide as a function of the electric field. Parameters used are / — 1500A, Bt = 0.03T, 
m = 17amu, A = 27r x 1.66GHz, /Xe = 1.667D, implying j/o — 10mm and b = 3T/m, and Y = lOOyo- 

meaning that the operation C simply changes the sign of the Hamihonian. A matrix representation for C can readily 
be found for the diagonal form of our molecular Hamiltonian 
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A physical interpretation of C remains to be found. 

D. Guide depth as a function of the electric field 

The guide depth can be found from Eqs. (5) and (7) of the main article which can be used to find the Stark-Zeeman 
energy E^/2,f{^Ty)- The guide depth can then be defined as 

Do = S3/2,/(0, Y) - ^3/2,/^, yo), (20) 

where F ^ j/o is a coordinate far away from the guide center. A plot of Do is shown in Fig[T]of this Supplement. As 
can be seen from the figure, for the field value at which anisotropy is exactly canceled [E ~ 3.5kV/cm, see Fig. 2(b) 
in the main article], the trap depth is larger than 2.7mK, which is the minimum OH temperature required for Eq.(2) 
of the main article to be valid. 




0\ 







1 

10 

1 

10/ 







(19) 



